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In metallic manganites with low Curie temperatures, a peculiar softening of the magnon spectrum 
close to the magnetic zone boundary has experimentally been observed. Here we present a theory of 
the renormalization of the magnetic excitation spectrum in colossal magnetoresistance compounds. 
The theory is based on the modulation of magnetic exchange bonds by the orbital degree of freedom 
of double-degenerate Cg electrons. The model considered is an orbitally degenerate double-exchange 
system coupled to Jahn- Teller active phonons which we treat in the limit of strong onsite repulsions. 
Charge and coupled orbital-lattice fluctuations are identified as the main origin of the urmsual 
softening of the magnetic spectrum. 
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I. INTRODUCTION 

The motion of charge carriers in the metallic phase 
of manganites establishes a ferromagnetic interaction be- 
tween spins on neighboring sites. Acppi'ding to the con- 
ventional theory of double exchangeJJTa the spin dynam- 
ics of the ferromagnetic state that evolves at temper- 
atm'es below the Curie temperature Tc is expected to 
be of nearest-neighbor Heisenberg type. This picture 
seems to be indeed reasonably accurate for manganese 
oxides with large values of Tq, i.e., for compounds whose 
ferromagnetiCpinetallic phase sustains up to rather high 
temperatures.Q However, recent experimental studies in- 
dicate marked deviations from this canonical behavior in 
compounds with low values of Tc- Quite prominent in 
this respect are measurements of the spin dynamics 
the ferromagnetic manganese oxide Pro.esSro.sTMnOa 
While exhibiting conventional Heisenberg behavior at 
small momenta, the dispersion of magnetic excitations 
(magnons) shows curious softening at the boundary of 
the Brillouin zone. This observation is of high impor- 
tance as it indicates that some specific feature of mag- 
netism in manganites has yet to be identified. 

A comparison of the magnetic behavior of different 
manganese oxides further highlights the shortcomings of 
the double-exchange theory: Assuming the magnon dis- 
persion to be of Heisenberg type, a small-momentum fit 
to a quadratic dispersion relation tOq — Dq^ yields the 
spin- wave stiffness D; in a conventional Heisenberg sys- 
tem the spin-wave stiffness scales with the strength of 
magnetic exchange bonds D oa J. Since the latter also 
controls the Curie temperature Tc ex: J, the ratio of D 
and Tc is expected to be a universal constant. Mangan- 
ites, on the other hand, exhibit a pronounced deviation 
from this behavior: D /Tc increases significantly as one 
goes from compounds with high to compounds with low 
values of Tea The presence of an additional mechanism 
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FIG. 1. The Bg-electron transfer amplitude, which controls 
the double-exchange interaction Jde, strongly depends on the 
orbital orientation: Along the z direction, e.g., d^22_^2 elec- 
trons (left) can hop into empty sites denoted by a sphere, 
while the transfer of d^2_„2 electrons (right) is forbidden. 



that controls the magnetic behavior of manganites is to 
be inferred. 

In the present paper, we propose a mechanism to ex- 
plain the above peculiar magnetic properties of ferromag- 
netic manganites. Our basic idea is the following: The 
strength of the ferromagnetic interaction at a given bond 
strongly depends on the orbital quantum number of Cg 
electrons (see Fig. ||) — along the z direction, e.g., only 
electrons in d'i,^2_^i orbitals can hop between sites and 
hence can participate in double-exchange processes; the 
transfer of (1^2 _y2 electrons is blocked due to the van- 
ishing overlap with O 2p orbitals located in-between two 
neighboring Mn sites. Temporal fluctuations of Cg or- 
bitals may thus modulate the magnetic exchange bonds 
(see Fig. 0), thereby renormalizing the magnon disper- 
sion. Short- wavelength magnons are most sensitive to 
these local fluctuations and are affected most strongly. 
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FIG. 2. Fluctuation of magnetic exchange bonds: Full lines 
denote active bonds, dashed lines inactive ones. 



Quantitatively the modulation of exchange bonds is con- 
trolled by the characteristic time scale of orbital fluctu- 
ations: If the typical frequency of orbital fluctuations is 
higher than the one of spins fluctuations, the magnon 
spectrum remains mostly unrenormalized — the orbital 
state then effectively enters the spin dynamics only on 
time average which restores the cubic symmetry of ex- 
change bonds. On the other hand, if orbitals fluctuate 
slower than spins, the renormalization of the magnon 
spectrum is most pronounced — the anisotropy imposed 
upon the magnetic exchange bonds by the orbital de- 
gree of freedom now comes into play. The presence of 
Jahn- Teller phonons enhances this effect by quenching 
the dynamics of orbitals. The suppression of fluctua- 
tions becomes almost complete as orbitals begin to or- 
der, resulting in a distinct softening of magnons which 
we interprete as a precursor effect of static orbital order. 
In the following, we calculated the dispersion of one- 
magnon excitations at zero temperature. We start from 
an orbitally degenerate Hubbard model that comprises 
the strongly correlated nature of the Mn 3d electrons 
and the physics of double exchange. The metallic motion 
of charge carriers establishes magnetic double-exchange 
bonds which are found to be further contributed to by 
virtual superexchange processes. Both types of exchange 
interaction are of ferromagnetic nature in the orbitally 
degenerate system subject to a strong Hund's coupling. 
Employipg a 1/5" expansion of spin and an orbital-liquid 
schemeDiiy to handle correlation effects, three different 
mechanisms are analyzed with respect to their capabil- 
ity to renormalize the magnon spectrum: scattering of 
magnons on orbital fluctuations, on charge fluctuations, 
and on phonons. Within this picture we can successfully 
reproduce the experimentally observed softening of the 
magnon dispersion. Furthermore we predict the renor- 
malization effect to become dramatic as static order in 
the orbital-lattice sector is approached. We note that 
the renormalization of the magnetic excitation spectrum 
by opticaLjjhonons has recently been investigated by 
Furukawa.EiJ 



II. MAGNETIC EXCHANGE BONDS 



The main aspects of the physics of manganites, i.e., 
the correlated motion of itinerant Cg electrons and the 
ferromagnetic interaction of Cg spins with a background 
of localized core spins, is captured by the following or- 
bitally degenerate Hubbard model: 
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with P — {siaSifj + j). The first term in Eq. (nl) describes 
the intersite transfer of electrons within degenerate eg 
levels. Here, c|^^ creates an eg electrons with spin and 
orbital quantum numbers s and a//3, respectively. The 
spatial direction of bonds is specified by 7 S {x,y,z\. 
One of the important features of the orbitally degener- 
ate mod|el| is the nondiagonal structure of the transfer 
matricesta 
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where we have chosen a representation with respect to the 
orbital basis {jSz^ — r^), \x^ — y^)}. The second term of 
Eq. (||) describes the Hund's coupling between the itin- 
erant Cg electrons and the localized core spins S'^; the 
magnitude of this coupling is Jh- The spin operator s^a 
acts on eg electrons in orbitals a, while s, = ^^^ Sia de- 
notes the total eg spin at a given site. Finally, the last two 
terms in model (0) account for the intra- (inter-) orbital 
Coulomb interaction U ([/') and the Hund's coupling be- 
tween eg electrons in doubly occupied states, riisa is the 
number operators of eg electrons in the state defined by 
s and a, and nta = X^s "■*"■ Double counting is excluded 
from the primed sum in the last term of Eq. ([l|) . 

In analogy to the transformation from a conventional 
Hubbard to t-J model, Eq. (||) can be projected onto 
the part of the Hilbert space with no double occupan- 
cies in the limit of strong onsite repulsions U ^ t and 
{If — Jh) 3> i. Doubly occupied states are then allowed 
only in virtual superexchange processes. Due to the pres- 
ence of Hund's coupling, the energy level of these virtual 
states depends on the spin orientation of_core and eg 
spins — a rich multiplet structure follows £3 The prob- 
lem considerably simplifies in the limit of large Hund's 
coupling U, U' 3> Jh 3> t which we believe to be realistic 
to manganites: Transitions to the lowest-lying interme- 
diate state with energy Ui — U' — Jh in which core and 
eg spins are in a high-spin configuration then dominate; 
doubly occupied sites with different spin structures lie 
higher by an energy of the order of ex Jh and can be ne- 
glected. We hence obtain the following t-J Hamiltonian: 
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The first two terms in Eq. (g) describe the double- 
exchange mechanism in the limit of strong onsite repul- 
sions. All double occupancies of Cg electrons are pro- 



jected out by the constrained operators c. 
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which act only on empty sites. The third term in Eq. (H) 
describes the superexchange interaction between singly 
occupied sites. The strength of this interaction is con- 
trolled by JsE = (2tV[/i)[S'(2S'-(- l)]-\ where S denotes 
the total onsite spin of Zd electrons. It is important to 
note that in the present model with large Jh, superex- 
change is of ferromagnetic nature. This stems from the 
fact that Hund's coupling forbids any double occupancy 
of a single Cg orbital. Pauli's exclusion principle, which 
is responsible for the antiferromagnetic nature of con- 
ventional superexchange, is therefore ineffective in dic- 
tating the spin structure of the virtual state. Rather, the 
spin orientation in the intermediate state is controlled by 
Hund's coupling which favors a ferromagnetic alignment 
of spins. The superexchange term in Eq. (@) exhibits 
yet another peculiar feature: The amplitude of superex- 
change processes depends on the orbital states of the eg 
electrons involved. This information enters via the or- 
bital pseudospin operators 
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where the Pauli matrices a^ act on the orbital sub- 
space; the factor (-j — t^t^) in Eq. (g) accounts for the 
specific nondiagonal structure of the transfer matrices 
t"^ and ensures that no double occupancy of a single or- 
bital occurs which would be forbidden by Pauli's exclu- 
sion principle and the large Hund's coupling. We finally 
note that superexchange processes in an orbitally degep-i 
erate system have also been studied by Feiner and Oles.E^ 
In the limit of large Jh the expression obtained by these 
authors maps onto the superexchange term of Eq. (0) for 
the special case S = 2. 

In the following, double-exchange and superexchange 
interactions which are jointly responsible for ferromag- 
netism in manganites are discussed in more detail. 



A. Double-Exchange Bonds 



We begin by analyzing the kinetic term of Hamiltonian 
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FIG. 3. The itinerant Cg spin (top left) interacts with the 
localized core spins (bottom left) via Hund's coupling. In the 
limit Jh S> t, the former can be separated from the orbital 
and charge degrees of freedom of the eg electron (circle) and 
can be absorbed into the total spin (bottom right). 



which establishes the double-exchange mechanism in the 
correlated system. Due to the strong Hund's coupling, 
core spins S'^ and itinerant eg spins s are not independent 
of each other; rather a high-spin state with total onsite 
spin S — S'^ + ^ is formed. This unification of band and 
local spin subspaces suggests to decompose the eg elec- 
tron into its spin and orbital/charge components. The 
eg spin can then be absorbed into the total spin, allow- 
ing an independent treatment of spin and orbital/charge 
degrees of freedom (see Fig. ||). The procedure of this 
separation scheme is the following: In a first step we in- 
troduce Schwinger bosons d,j| and dii (see, e.g., Ref. |lj) 
to describe the Cg spin 
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as well as Schwinger bosons dL and dJ, to model the 
total onsite spin 



s+ = dId 



sr = dI d, 



These auxiliary particles are subject to the following con- 
straints that depend on the Cg occupation number n^: 
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The creation and destruction operators for eg electrons 
can then be expressed in terms of spinless fermions Cia 
which carry charge and orbital pseudospin and Schwinger 
bosons which carry spin: 



^t = - E E<' (4aS,s;3 + H.C.) - Jh E ^.'^- 
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The kinetic-energy Hamiltonian (^ now describes the 
transfer of pairs of spinless fermions and Schwinger 
bosons: 
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The Bose operators are subject to the constraint (0) that 
enforces the operators dis and dj^ to act only on pro- 
jected Hilbert spaces with one or zero Schwinger bosons, 
respectively. Our aim is to absorb the eg spin into the 
total spin, which requires to map the eg operators dis 
onto operators Dis for the total spin. This is done by 
comparing the matrix elements of the two types of op- 
erators. On the one hand, keeping in mind that Hund's 
rule enforces the onsite spins to be always in a total-spin- 
symmetric state, the only nonvanishing matrix elements 
of the dis operators are 
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S,m) = ^iS + m)/{2S), (7) 

S,m) = ^iS-m)/i2S). (8) 



In deriving the above expressions we have used the 
Clebsch-Gordan coefficients {S'^ ,rrf]rrf\S^m) to decom- 
pose the total-spin state \S,m) into core- and eg-spin 
states |S"^,m'^;m^) with mf' =| / [. These coefficients 
are given by 
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On the other hand, the matrix elements of the Dis oper- 
ators are 
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S, m^ ^ ^{S + m), (9) 

S,mj = ^/{S-m). (10) 



All other matrix elements vanish due to the constraint 
of Eq. d). By comparing Eqs. (0)-(|) with Eqs. (|)-(|lC 
we obtain the mapping 
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Hamiltonian (|6|) can hence be rewritten in terms of total- 
spin operators Dis'. 
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The Hund's coupling term of Eq. (o) has been dropped 
here as its presence is implied by the spin construction 
employed above. This completes the separation of spin 
from the charge/orbital quantum numbers of eg elec- 
trons. 



At low temperatures the magnetic moment of ferro- 
magnetic manganites studied here is almost fully sat- 
urated. It is therefore reasonable to expand Eq. ( [ll| ) 
around a ferromagnetic groundstate. Technically this is 
done by condensing the spin-up Schwinger bosons (as- 
suming the ferromagnetic moment to point along this 
direction) and by treating spin-wave excitations around 
this groundstate in leading order of 1/S. Introducing 
magnon operators 6^, the following relations hold: 
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This spin representation fixes the number of Schwinger 
bosons per site to 25. The essence of the 1/5 expansion is 
to consider the presence of a hole as a small perturbation 
which changes the spin projection 5^ but not the spin 
magnitude 5. Employing magnon operators, the kinetic- 
energy Hamiltonian (pT|) hence becomes 
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The first term of Eq. ( [l^ ) describes the motion of strongly 
correlated fermions in a ferromagnetic background. The 
second term controls the dynamics of spin excitations 
in the magnetic background and the interaction of these 
excitations with the fermionic sector. 

At small magnon numbers, i.e., at low temperatures 
T <C Tc, Eq. (|lj) can be mapped onto the following 
expression for the magnetic double-exchange bonds: 
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Equation ( |13D highlights an important point: The 
strength of double-exchange bonds is a fluctuating com- 
plex quantity. Only when treating the orbital and charge 
sectors on average, i.e., when replacing the bond opera- 
tors cj^Cip by their mean-field value {Ci^Ci/3), an effective 
Heisenberg model as in a conventional mean-treatment of 
double exchange is obtained: H = Jde J^uj) '^i'^j with 

Jde = (252)-i J2c.f3^ tficLcjp)- In Section || we inves- 
tigate in more detail the modification of the mean-field 
picture by fluctuations in the bond amplitude. 

It is interesting to turn to the limit of classical spins 
shortly. Replacing the spin operators in Eq. (|l3|) by their 
classical counterparts 5^ = 5 cos 6* and 5* — 5 sin 6*6^*"^, 
an effective fermionic model is obtained: 
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This model exhibits-an unconventional phase-dependent 
hopping aniplitude:cJ 
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A similar result has been discussed in Refs. |16|,|17|in terms 
of a Berry-phase effect. 



B. Superexchange Bonds 

At low- and intermediate-doping levels, virtual charge- 
transfer processes across the Hubbard gap becomes of 
importance. These superexchange processes establish an 
intersite interaction, which in the limit of a strong Hund's 
coupHng is described by [see Eq. (ph] : 
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As mentioned above, superexchange is of ferromagnetic 
nature in the orbitally degenerate system with strong on- 
site correlations. Double exchange and superexchange 
therefore act together in establishing |tijc ferromagnetic 
exchange links in metallic manganites.ta 

Following the discussion on double-exchange bonds we 
express the spin operators in Eq. dl3) in terms of magnon 
operators bi. This leads to 
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Equation (n6h describes the interaction between orbital 
fluctuations and the magnetic sector of the Hilbert 
spacclij The phase dependence exhibited by the double- 
exchange counterpart Eq. (^ is absent here. This is due 
to the fact that superexchange is a second-order process 
which depends only the amplitude but not on the phase 
of the transfer amplitude. 



III. MAGNON DISPERSION 

In the previous section, the role of double-exchange 
and superexchange processes in promoting ferromag- 
netic exchange bonds in manganites was discussed. At 
intermediate-doping levels these exchange interactions 
induce a ferromagnetic groundstate in a variety of man- 
ganese oxides. We now turn to analyze the propagation of 
magnetic excitations in this ferromagnetic phase, namely 
by deducing the dispersion relation of single-magnon ex- 
citations. 



In a first step, we derive the correct operator for creat- 
ing a magnetic excitation in hole-doped double-exchange 
systems. It has to account for the fact that the total on- 
site spin depends on whether a hole or an eg electron is 
present at that site: The spin number is S* — 2 in the for- 
mer and S in the latter case. This difference in the spin 
number was neglected in the 1/S expansion employed 
in Sec. ^ Here this approximation is no longer valid, 
which requires a rescaling of the magnon operators bi. In 
general, a spin excitation is created by the operator S^ . 
Expressing this operator in terms of Schwinger bosons 
S^ = DJ^Dii, condensing Di^, and mapping Di^ onto 
the magnon operator bi, the following representation is 
obtained: 



S+ = 



2S bi, for sites with Cg electron, 

VSS* — 1 bi, for sites with hole. 



Assuming S to be the "natural" spin number of the sys- 
tem, the magnon operator bi hence has to be rescaled by 
a factor [(25" — 1)/(2S')]^/^ when being applied to hole 
sites: 
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for sites with e„ electron, 
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The general magnon operator that automatically probes 
the presence of an Cg electron can finally be written as 
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where rii is the number operator of Cg electrons. Bi repre- 
sents the true Goldstone operator of hole-doped double- 
exchange systems. Its composite character comprises lo- 
cal and itinerant spin features which is a consequence of 
the fact that static core and mobile eg electrons together 
form the total onsite spin. While the itinerant part of 
Bi is of order 1/5 only, it nevertheless is of crucial im- 
portance to ensure consistency of the spin dynamics with 
the Goldstone theorem, i.e., to yield an excitation mode 
whose energy vanishes at zero-momentum. 

Having derived the correct magnon operator for doped 
double-exchange systems, we now study the propagation 
of the magnetic excitations it creates. The link between 
sites that allows a local excitation to spread throughout 
the system is established by the exchange-bond Hamilto- 
nians (O) and dlo) . At low temperatures the dynamics of 
spin waves which hence develop is captured by the single- 
magnon dispersion. The important question we are inter- 
ested in is the following: To which extent is the magnon 
spectrum affected by fluctuations in the exchange bonds? 
To answer this question we express the full magnon spec- 
trum (Dp in terms of the conventional mean-field disper- 
sion LOp and the magnon selfenergy Y,(lj,p): 
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Fluctuation are considered only on average in the former 
but are explicitly accounted for in the latter term. The 
mean-field dispersion Wp as well as the scattering vertices 
needed to construct S(w, p) can be derived by commuting 
the magnon operator Bi with the Hamiltonian. To be 
specific we explicitly perform this commutation, for now 
restricting ourselves to the double-exchange Hamiltonian 
Ht given by Eq. (|l2|). In the momentum representation 
we obtain 

[Bp, Ht] — ijJpBp 
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The two terms on the r.h.s. of Eq. (O) correspond to 
an expansion of the bond operators cl^Cjp around their 
average value: 
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The mean-field magnon dispersion tOp in the first term of 
Eq. (|8|) is of conventional nearest-neighbor Heisenberg 
form 
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with the form factor 7p = z~^ ^^ exp(ip^), z — 6, and 
the spin- wave stiffness constant is D = SJue- On this 
mean-field level the strength of the exchange bonds de- 
pends on the orbital and charge degrees of freedom only 
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ond term in Eq. (|18| ) is the scattering vertex needed to 
construct the magnon selfenergy I](ijj, p). It describes the 
interaction between magnons and orbital/charge fiuctu- 
ations. The vertex function is 



Af{k) = 7f - 7.t, 



with the form factor 7^ — (zt) ^ ^^ 
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vertex function Ap^{k) vanishes in the limit p 
compliance with the Goldstone theorem. 

Before we can engage in evaluating the magnon selfen- 
ergy associated with the scattering vertex in Eq. (nS) , the 
problem of dealing with the correlated nature of fermionic 
operators cj^ = ^^(1 — rii) has to be addressed. To han- 
dle the constraint that allows only for one, e_ electron per 
site, we employ an orbital-liquid scheme:BM Orbital and 
charge degrees of freedom of the Cg electron are treated on 
separate footings by introducing "orbiton" and "holon" 
quasiparticles. To describe an orbitally disordered state 
in which orbitals fluctuate strongly, orbitons fi axe as- 
signed fermionic and holons hi bosonic statistics.EEl The 
original fermion operators are hence reexpressed by 
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The local no-double-occupancy constraint is now relaxed 
to a global one: 
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The main feature associated with the constrained nature 
of electrons, namely the separation of energy scales of 
orbital and charge dynamics, sustains this procedure due 
to the fact that two different types of quasiparticles are 
being used. Introducing mean-field parameters 
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where x is the concentration of holes in the system, or- 
bitons and holons can now be decoupled. We note that 
the two mean-field parameters in Eq. ( pi[ ) are approxi- 
mately related by x = ^(1 — xY 

Employing representation (EG), we reexpress the com- 
mutator of Eq. ( |lq ) in terms of orbiton and holon oper- 
ators: 
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The vertex functions are given by 
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Orbitons and holons have been decoupled in Eq. (p2) 
by employing the mean-field parameters x and x of Eq. 
(pl|). This yields two different types of scattering vertices, 
one describing the interaction of magnons with orbital 
ffuctuations, i.e., orbitons, the other of magnons with 
charge fluctuations, i.e., holons. 

Finally we include in our treatment the magnetic 
bonds stemming from superexchange processes as de- 
scribed by Hj in Eq. (p^. The effect is twofold: Su- 
perexchange enhances the spin-wave stiffness D which 
now becomes 

D - 5(Jde + Jse) = tx{x + a;o)/(25), 

with xq — 2xt/Ui; further, superexchange processes 
renormalize the vertex function of magnon-orbiton scat- 
tering which becomes 

Cf{k)^xAf{k)+xoBf{k,q), 



with 



_ Q/3 



Q/3 



Q/3 



Q/3 



5p'^(fc, g) = 7r + Ik'-q - 7.Tp - <-q- 



From the two types of scattering vertices in Eq. ( |2^ ) , two 
contributions to the magnon selfenergy follow. These 
describe the scattering of magnons on orbitons and on 
holons and are depicted in Figs, ^a) and (b), respec- 
tively. 




(a) 



(b) 



(c) 



FIG. 4. Magnon selfenergies describing the effect of 
magnon scattering on (a) orbital fluctuations, (b) charge fluc- 
tuations, and (c) phonons. Solid, dashed, dotted, and wiggled 
lines denote orbiton, holon, magnon, and phonon propagators, 
respectively. 



An important piece of physics is still missed in the 
above treatment, namely the Jahn- Teller coupling of or- 
bitals to the lattice£3 In a cubic system there exist two 
independent Jahn- Teller modes Q2 and Q^ which lift the 
degeneracy of singly occupied Cg orbitals. The interac- 
tion between orbitals and these two orthogonal lattice 
modes is described by 



H 



JT 



E(52Q 



2i(^'i + gsQsicrt 



X j z 



(23) 



where the Pauli matrices a^ act on the orbital sub- 
space and the coupling constants 172 ~ 53- The crystal 
dynamics is controlled by the Hamiltonian 
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ph 
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(24) 



l^j/y 



with Q'l'" = {Q3, ± VSQ2^)/2, Ql = Q^,, and Q, = 
(Q2i,Q3i); Pi denotes the conjugate momentum vector 
corresponding to the lattice distortions Qi. The three 
terms on the r.h.s. of Eq. (e3) account for the crystal 
deformation energy, the correlations between neighbor- 
ing sites, and the lattice kinetics, respectively. Equation 
( P4| ) can be diagonalized in the momentum representa- 
tion, yielding 



ffph - / ^^kO-ku'^kv, 



(25) 



hu 



with index i/ = ± and the phonon dispersions 

1/2 



uj^ =0)0 Ktfc ± y n 



''2k ^ "-Sfc 



(26) 



A 1 (2) 

Cz), K2fe = Kl% , K3fc 

= -V3(Ca; - Cy 



)/2, 



^Cy with Cq — cos ka, and wq — ^J K/M. 



Here, Kifc = 1 + ki{cx + Cy 

^iVk with fci = Ki/K and 77^" 

(3) r 

% -C.-2 

While there is no direct coupling between spins and 
phonons in the present system, lattice modes never- 
theless strongly affects the spin dynamics. The link 
between spin and lattice is established via the orbital 
channel: The coupling of orbitals to the lattice imposes 
low phononic frequencies onto orbital fluctuations. This 
acts to enhance the modulation of magnetic exchange 





FIG. 5. Effective spin-phonon-coupling vertex. The domi- 
nant contribution shown on the right stems from a combina- 
tion of spin-orbital- (filled dot oc t) and orbital-lattice- (open 
dot oc g2) coupling vertices. The orbital susceptibility de- 
picted by a bubble controls the coupling strength. Solid, dot- 
ted, and wiggled lines represent orbiton, magnon, and phonon 
propagators, respectively. 



bonds; thereby the effect of phonons extends onto the 
spin sector. To study this mechanism in more detail, we 
construct an effective spin-phonon-coupling Hamiltonian 
from which we then calculate the phononic contribution 
to the magnon selfenergy. Combining the spin-orbital- 
coupling term of the exchange Hamiltonians (03) and 
( [l6|) with the orbital-lattice Hamiltonian (|2^) we obtain 
(see Fig. |): 



Hs-ph - - E 9pqi%i, + 0.qu)BlBpj^ 
pqu 

The coupling constants in Eq. (E^) are 

1/2 
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= £0 — Al?) sin e„ + All) cos 6 



pq 



pq 



1 ' 



with eo = {E,Talcoo/S^y/' and A^f,' = (r;'") - 4"^ 
Vp+q)- Further 
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v2 



sineq = -= I 1 



'*2q + '^3q 



K3q 
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sign(K2q). 
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The strength of the spin-lattice interaction is controlled 
by the orbital susceptibility {{fj^z tfi'\)i'^i))(^ which en- 
ters the parameter oq = t{x+xo){ifi^^^-^ fi^){al))uj=o; ^^^ 
zero-frequency limit is admissible bearing in mind that 
the energy scale of orbital fluctuations exceeds the one 
of phonons. The phononic contribution to the magnon 
selfenergy that follows from Hamiltonian (|2^) can finally 
be calculated, the corresponding diagram is depicted in 
Fig. 1(c). 
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FIG. 6. Magnon dispersion along (0,0,^), {(,,(,,0), and 
{£,, ^, ^) directions, where ^ = 0.5 at the cubic zone bound- 
ary. Experimental data from Ref. in are indicated by circles, 
the mean-field dispersion tOp of Eq. ^19) is marked by a dashed 
line; the latter is of conventional nearest-neighbor Heisenberg 
form. Solid lines represent the theoretical result for the dis- 
persion ujp defined by Eq. (|l7]); it includes charge, orbital, 
and lattice effects. The upper curve is obtained for disper- 
sionless phonons with fei = 0, the lower one is a fit to the ex- 
perimental data with fci = —0.33 corresponding to ferrotype 
orbital-lattice correlations. 



IV. COMPARISON WITH EXPERIMENT 

We are now in the position to evaluate the selfenergies 
of Fig. ^. Charge and orbital susceptibilities are calcu- 
lated using mean-field Green's functions in slave-boson 
hi and fermion fi subspaces. For the spectral density of 
Jahn- Teller phonons in Fig. H(c) we employ the expres- 
sion 



.Ph 



1 



P± (^, 9) = - 



7rc^±(a;-cj±)2+r2^ 



(28) 



which phenomenologically accounts for the damping F 
of phonons due to their coupling to orbital fluctuations. 
The phonon dispersion lu^ is given by Eq. (|2q). 

The expressions obtained from the diagrams in Fig. 
y contain summations over momentum space which we 
perform numerically using a Monte-Carlo algorithm. The 
result is shown by solid lines in Fig. |6|. For comparison, 
the experimental data of Ref. M are marked by circles and 
the bare mean-field dispersion ujp is indicated by a dashed 
line. The following parameters are chosen: The hopping 
amplitude t = 0.4 e\Lis adjusted to fit the spin stiffness 
in Pro.esSro.sTMnOaP further we use Ui = 4 eVJB The 
phonon contribution depends on the quantities Ejt^Uq = 



{g2aof/2K ^ 0.004 eVj^ wq = 0.08 eVP and F 
0.04 eV. 

The upper solid line in Fig. g is obtained for fci = 0. In 
this case intersite orbital-lattice correlations in Hamilto- 
nian (e3) are discarded — phonons are dispersionless. A 
pronounced softening of magnons at large momenta can 
be observed. A more detailed analysis reveals this effect 
to be mostly due to fluctuations of the orbital and lattice 
degrees of freedom. In contrast, charge fluctuations are 
found to play only a minor role. We attribute this to the 
fact that the spectral density of charge fluctuations lies 
well above the magnon band. Orbital and lattice fluc- 
tuations, on the other hand, are of rather low frequency 
(ex xt and oc Wq , respectively) and hence affect the spin- 
wave dispersion in a more pronounced way. 

The lower solid line in Fig. ^is obtained for fci = —0.33 
which yields a fit to the experimental data of Ref. [j. 
The directional dependence of the magnon renormaliza- 
tion seen in experiment is well reproduced: The effect 
is strongest in (0,0,^) and (0,^,^) directions. A key 
observation here is the crucial role of intersite correla- 
tions of orbital-lattice distortions — these are captured 
by the phononic dispersion being controlled by the pa- 
rameter fci . In order to reproduce the experimental data 
we are forced to assume these correlations to be of fer- 
rotype^ i.e., fci < 0. We believe this somewhat surprising 
result to reflect an important piece of new physics: Con- 
ventionally one would expect fci > associated with a 
tendency of the orbital/lattice sector to develop antifer- 
rotype orderJij In the hole-doped system, however, this 
effect competes against charge mobility which prefers a 
ferrotype orbital orientation. The latter allows to mini- 
mize the kinetic energy by maximizing the transfer am- 
plitude between sites. While Jahn- Teller lattice effects 
prevail at low doping, we speculate the kinetic energy to 
dominate at large enough hole concentrations. In fact, 
low-dimensional ferrotype orbital correlations (resonat- 
ing \x^ — j/2), \x'^ — z^), \'ip — z^) planar configurations) 
have been observed to, evolve in a bosonic description 
of orbital fluctuations.cl The fermionic description of or- 
bitals employed in the present work emphasizes on mod- 
eling a strongly fluctuating orbital-liquid state, but un- 
derestimates these orbital-lattice instabilities. In order 
to simulate the competition between Jahn- Teller effect 
and kinetic energy we therefore turn to a phenomenolog- 
ical approach: By tuning the parameter fci we control 
the character of intersite orbital-lattice correlations. The 
result for different values of fci is shown in Fig. M where 
SjTog = 0.006 eV is used. 

Ferrotype orbital correlations with fci < arc found to 
be most effective in renormalizing the magnon spectrum. 
This is ascribed to slowly fluctuating layered orbital con- 
figurations which effectively reduce the dimensionality of 
exchange bonds. We note that magnons in (^, ^, ^) di- 
rection are sensible to all three spatial directions of the 
exchange bonds; their dispersion therefore remains unaf- 
fected by the local symmetry breaking induced by low- 
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FIG. 7. Magnon dispersion including charge, orbital, and 
lattice effects. Different values for fci controlling intersite or- 
bital-lattice correlations are used. The softening enhances as 
fci — > — ^ corresponding to an instability point towards fer- 
rotype orbital-lattice order. 



dimensional orbital correlations. As an instability to- 
wards orbital-lattice order is approached, exchange-bond 
fluctuations become quasistatic. In the magnon spec- 
trum this is reflected by a strong enhancement of the 
renormalization effect as is seen in Fig. for fci -^ — i. 
The layered orbital structure which evolves at this point 
is accompanied by a layered spin structure; the latter is 
indeed etnerimentally observed at doping levels of about 

We finally note that the softening of magnons at the 
zone boundary leads to a reduction of Tc- Remarkably, 
the small- g spin stiffness D remains unaffected which ex- 
plains the anomaloiis enhancement of the D /Tc ratio in 
low- Tc manganites-Q 



V. CONCLUSION 

In summary, we have presented a theory of the spin 
dynamics in ferromagnetic manganites. Taking into ac- 
count the orbital degeneracy and the correlated nature 
of Cg electrons, we analyzed the structure of magnetic 
exchange bonds; these are established by the intersite 
transfer of electrons in coherent double-exchange and vir- 
tual superexchange processes. Orbital and charge fluc- 
tuations are shown to strongly modulate the exchange 
bonds, leading to a softening of the magnon excitation 
spectrum close to the Brillouin zone boundary. The pres- 
ence of Jahn- Teller phonons further enhances the effect. 
This peculiar interplay between double-exchange physics 
and orbital-lattice dynamics becomes dominant close to 
the instability towards an orbital-lattice ordered state. 
The unusual magnon dispersion experimentally observed 
in low- Tc manganites can hence be understood as a pre- 
cursor effect of orbital-lattice ordering. While the soft- 
ening of magnons at the zone boundary is responsible for 
reducing the value of Tc , the small-momentum spin dy- 



namics that enters the spin- wave stiffness D remains vir- 
tually unaffected. This explains the enhancement of the 
ratio T) jTc observed in low- Tc compounds. In general 
it can be concluded that strong correlations and orbital 
fluctuations play a crucial role in explaining the peculiar 
magnetic properties of metallic manganites. 
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